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O ■ ABSTRACT 

^ ' Elucidating the processes that governed the assembly and evolution of galaxies over cosmic 

^ , time is one of the main objectives of all of the proposed Extremely Large Telescopes (ELT). 

To make a leap forward in our understanding of these processes, an ELT will want to take ad- 
vantage of Multi-Objects Adaptive Optics (MOAO) systems, which can substantially improve 
the natural seeing over a wide field of view. We have developed an end-to-end simulation to 
^ ' specify the science requirements of a MOAO-fed integral field spectrograph on either an 8m 

! or 42m telescope. Our simulations re-scales observations of local galaxies or results from 

On . numerical simulations of disk or interacting galaxies. The code is flexible in that it allows 

us to explore a wide range of instrumental parameters such as encircled energy (EE), pixel 
size, spectral resolution, etc. For the current analysis, we limit ourselves to a local disk galaxy 
OO ' which exhibits simple rotation and a simulation of a merger. While the number of simulations 

' is limited, we have attempted to generalize our results by introducing the simple concepts of 

OO . "PSF contrast" which is the amount of light polluting adjacent spectra which we find drives 

■ the smallest EE at a given spatial scale. The choice of the spatial sampling is driven by the 

"scale-coupling". By scale-coupling we mean the relationship between the IFU pixel scale and 
the size of the features that need to be recovered by 3D spectroscopy in order to understand 
/\ ' the nature of the galaxy and its substructure. Because the dynamical nature of galaxies are 

mostly reflected in their large-scale motions, a relatively coarse spatial resolution is enough to 
distinguish between a rotating disk and a major merger Although we used a limited number 
of morpho-kinematic cases, our simulations suggest that, on a 42m telescope, the choice of 
an IFU pixel scale of 50-75 mas seems to be sufficient. Such a coarse sampling has the benefit 
of lowering the exposure time to reach a specific signal-to-noise as well as relaxing the per- 
formance of the MOAO system. On the other hand, recovering the full 2D-kinematics of z~4 
galaxies requires high signal-to-noise and at least an EE of 34% in 150 mas (2 pixels of 75 
mas). Finally, we carried out a similar study for a hypothetical galaxy/merger at z-1.6 with a 
MOAO-fed spectrograph for an 8m, and find that at least an EE of 30% at 0.25 arcsec spatial 
sampling is required to understand the nature of disks and mergers. 

Key words: Galaxies: evolution - Galaxies: high-redshift - Galaxies: kinematics and dynam- 
ics - Instrumentation: adaptive optics - Instrumentation: high angular resolution - Instrumen- 
tation: spectrographs 



1 INTRODUCTION 

Developing a coherent model for the mass assembly of galaxies 
over cosmic time is a complex and difficult task. It is difficult 
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because developing such a coherent model involves highly non- 
linear physics (e.g., the cooling and collapse of baryons, or feed- 
back from stars and super-massive blackholes which regulate both 
the growth of the stellar mass and black holes themselves) and all 
over a very wide range of physical scales (from large scale struc- 
ture to star clusters to black holes). Our knowledge of galaxy for- 
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mation in situ relies mostly on observations of integrated quan- 
tities,which is insufficient to understand the details of the pro- 
cess of formation, e.g., the interplay between the baryonic and 
dark matter, or the complex physics of baryons such as merging, 
star formation, feedback, various instabilities and heating/cooling 
mechanisms, etc. At the heart of our developing understanding of 
galaxy evolution is the relative importance of secular evolution 
through (quasi-)adi abatic accretion of m ass, instabilities, and res- 
onances (e.g., Semelin & CombeslllOOSi), versus more violent evo- 
lution through merging (e.g., Hammer et al.ll200^ . as a function of 
look-back time. 

Because of the complexity of the processes involved, under- 
standing the mass assembly of galaxies over cosmic time requires 
constraining a wide range of phenomenology for large samples of 
objects. In this respect integral field spectroscopy is particularly 
well suited in that it allows us to map the spatially resolved phys- 
ical properties of galaxies (see, e.g., Puech et al. 2006b). Unfortu- 
nately, splitting the light into many channels in a integral field spec- 
trograph requires long integration times, even on large telescopes, 
to detect low surface brightness emission. Thus, conducting large 
statistical studies necessary for understanding the processes driv- 
ing galaxy evolution will demand a high multiplex capability to 
efficiently investigate sufficient numbers of objects. 

Using current facilities on 8m class telescopes, it is already 
possible to map the physical properties of distant galaxies. For ex- 
ample, FLAMES/GIRAFFE on the VLT offers the unique ability to 
observe 15 distant galaxies simultaneously. The first (complete) 3D 
sample of galaxies at moderate redshifts has revealed a large frac- 
tion (~40%) of all z~0.6 galaxies with Mg <-19.5 have per turbed 
or complex kinematics (Flores et al. 20061: pYang et aLlFlOOSh . This 
population of dynamically disturbed galaxies is most probably a re- 
sult of a high preva l ence o f mergers and merger remnants (see also 
IPuech et al.ll2006al . |2008|) . At higher redshift, integral field spec- 
troscopy of several objects has been obtained using single object in- 
tegral fieM^peclrographs avai lable on 8-10 meter class telescopes 
(e.g.. iForster-Schreiber et alJi2006h . However, many of these data 
sets were obtained with limited spatial resolution sampling galax- 
ies on scales of a few kpc, and the true dynamical nature of many 
sources sampled at these spatial and spectral resolutions remains 
uncertain (e.g.. Law et al. 2007). 

Improving and extending this approach to the general under- 
standing of the nature of faint galaxies and to galaxies in the early 
Universe will require a new generation of multi-object integral field 
spectrographs working in the near-IR, with increased sensitivity 
and better spatial resolution on even larger telescopes. Because of 
the complex interplay between spatial and spectral features and our 
limited knowledge of high redshift galaxies, it is helpful, perhaps 
necessary, to rely on numerical simulations for constraining the de- 
sign parameters of this new kind of instrument. For this purpose, 
we have developed software that simulate end-to-end the emission 
line characteristics of local galaxies and numerical simulations of 
galaxies to show how they would appear in the distant Universe. 
By changing various parameters like the resolution, pixel scale, 
and point spread function it is possible to constrain the instrumen- 
tal characteristics and performance against a set of galaxy charac- 
teristics (e.g., velocity field). This paper is the first in a series to 
investigate the scientific design requirements of possible future in- 
struments, telescope aperture and design, and site characteristics 
appropriate for constrain how galaxies grew and evolved with cos- 
mic time. In this paper in particular, we aim to present the details of 
our assumptions that go into the software, as well as a first attempt 
to constrain some of the high level specification for a hypotheti- 



cal but realistic multi-object integral field spectrograph assisted by 
MOAO system. Because 2D kinematics (i.e., velocity field and ve- 
locity dispersion maps) are one of the most obvious and easiest to 
currently simulate, they provide the most straight forward way of 
constraining the high level science design requirements for this type 
of instrument. We emphasize that the goal of this paper is not to 
study the detailed scientific capability of such instruments. Instead, 
the present paper aims at examining a few scientifically motivated 
cases to derive relevant specifications for such types of instruments 
with a capable MOAO system. These specifications will then be 
adopted as a baseline for studying their scientific capabilities sub- 
sequent papers. This paper is organized as follows: in Section 2, 
we detail the goals of this study and our methodology: in Sect. 3 
we present the new simulation pipeline and in Sect. 4 how MOAO 
PSFs are simulated: In Sect. 5 we describe the kinematic measure- 
ments: In Sect. 6 and 7, we present the simulations and their results, 
which are discussed in Sect. 8. In Sect. 9, we draw the conclusions 
of this study. 



2 SPECIFYING MOAO-FED SPECTROGRAPHS 

2.1 MOAO & 3D spectroscopy: EAGLE and FALCON 

Most of current dynamical studies of distant galaxies are severely 
hampered by their limited spatial resolution, which introduces an 
important uncertainty as to their dynamical nature as discussed in 
the Introduction. A MOAO-fed integral field spectrograph on the 
VLT and/or the E-ELT will improve substantially the spatial reso- 
lution compared to non-AO assisted observations. The concept of 
MOAO was first proposed in the FALCON de sign study, a multi ob - 
jects 3D spectrograph intended for the VLT jHammer et al.ll2002l) . 
Briefly, the main objectives of the FALCON concept was to obtain 
3D spectroscopy of distant emission line galaxies, up to z ~ 2. Do- 
ing this requires us to observe in the near infrared since all of the 
important optical emission lines from [OII]/13727A to Ho- are red- 
shift ed into the atmospheric windows from 0.8 to 2.5 yum (see, e.g., 
, Hammer et al .112004,) . Because the intrinsic size of distant galax- 
ies is small, AO is needed to spatially resolve emission lines and 
thus to keep the spectroscopic signal-to-noise ratio (SNR) high 
enough to detect small features. Correcting a large field of view 
(FoV) to observe many galaxies simultaneously is challenging and 
normal single natural guide star AO is not sufficient. MOAO is a 
method for obtaining excellent correction over a relatively large 
FoV (>few arcmin). The basis of this concept is to correct only re- 
gions where the targets lie, instead of correcting the entire FoV of 
the instrument. MOAO uses atmospheric tomography techniques 
( lTallon&Fovlll99(]| :l Tokovinin et al.ll200lh meaning several wave- 
front sensors (WFS) for each integral field unit (IFU) measure the 
off-axis wavefront coming from stars located within the FoV, and 
the on-axis wavefront from the galaxy is deduced from off-axis 
measurements and corrected by an AO system within the optical 
tr ain of each IFU. A de tailed description of FALCON can be found 
in lAssemat et al. li liooTh . The general concept of MOAO was subse- 
quently adapted to the future European Extremely Large Telescope 
(E-ELT) in the WFSPEC study (Moretto et al. 2006), which was a 
precursor to the current EAGLE concept. One of the most impor- 
tant goals in setting the science design requirements for EAGLE is 
mapping the physical properties of galaxies up to z~5, in order to 
better understand the processes that drive the evolution of galaxies. 
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UGC 


Type z 


DC) 


Ms 


inc(deg) 


5253 


SA(rs)ab 0.004306 


4.168 


-20.8 


40 261.8 



Table 1. Main photometric and dynamical properties of UGC5253, used 
as a typical rotatin g disk in the simulations: Morphological type (from 
iGarrido et al]|2002l) . redshift (from CDS), absolute B-band rnagnitude, in- 
clination, and maximal velocity of rotation (from lGarrido et alj2002h . 



2.2 Goals of this study 

As a first study, we want to help derive high-level design specifi- 
cations for the EAGLE and FALCON MOAO-fed instruments by 
addressing the following: 

• Explore the necessary order of magnitude of integration times 
that will be required to reliably derive the kinematics of distant 
emission line galaxies, allowing their classification; 

• By selecting a simple rotating disk with little substructure and 
a merger simulation with two overlapping disks, explore the most 
relaxed requirements for the design of the instrument in terms of 
the necessary adaptive optics performance and coarseness of the 
pixel scale. To this end we have limited our simulations to spatial 
sampling of 50-75 mas pixel"' ; 

• Help constrain the requirements on the MOAO system by 
specifying its spectroscopic coupling (which is defined here as the 
Ensquared Energy [EE] per element of the spatial resolution, see 
Sect. 4). 

A detailed study of the scientific capabilities of these instru- 
ments is beyond the scope of this paper. Instead, we want to use 
a limited number of comparisons between scientifically motivated 
cases in order to help derive the most relaxed high-level specifica- 
tions. 



2.3 Methodology 

Flux is the zero-order moment of an emission line, while the ve- 
locity and width are the first and second moment, and higher order 
moments always have more relative uncertainty. For this reason, 
kinematics will set stringent requirements on the SNR necessary to 
fully characterize faint extended emission line gas. 

The characteristics of these moments depend on the underly- 
ing source of an emission line's excitation and kinematics. Because 
the most easily reached goal of MOAO-fed spectrographs will be 
to understand the kinematics of the emission line gas in galaxies, a 
natural kinematic reference case for such a study is a simple spiral 
rotating disk. In the following we will use Fabry-Perot observations 
of the lo cal rotating disk U GC5253 taken as part of the GHASP 
survey ( AmrametalJl2002l: see its main photometric and dynami- 
cal properties in Tab.[T) see also Fig.[T]( as such a reference case. 

Several different mechanisms have been proposed to explain 
how galaxies assembled their mass: major mergers (merging be- 
tween galaxies having similar masses), minor mergers (merging be- 
tween galaxies with a mass ratio larger than 5:1), and gas accretion 
(of cold or hot gas from the IGM or dark matter halo). For each of 
these processes, there are a large number of possible configurations 
or characteristics that need to be modeled or analyzed (e.g., mass 
ratio of merging galaxies, the relative angular momentum vectors 
of each merging galaxy and of the orbit, the temperature and den- 
sity distribution of the gas as it is being accreted, etc.). However, 
our goal here is not to conduct a systematic study of the scientific 



capabilities of MOAO-fed spectrograph in discriminating between 
all possible evolutionary mechanisms. Instead, we want to inves- 
tigate the surface brightness limits of such instruments and focus 
on the possibility to use a relatively coarse spatial resolution to re- 
cover large scale motions, whereby the process underlying the dy- 
namics may be divined (following for example iFlores et al]|2006l : 
IPuech et alj2006al:lYang et al.l2008b . Because we only want to help 
determine plausible design requirements relative to this goal, it is 
sufficient to explore the most extreme large-scale kinematic pat- 
terns, which are produced by major mergers. Indeed, during such 
processes, the disks of the progenitors are usually completely de- 
stroyed, which leads to complex kinematic patterns with abrupt 
changes in the velocity gradient(s). In the following, we will use a 
major merger extracted from an hydrodynami c simulation of a pair 
of two merging Sbc galaxies jCox et alj2006h . just before the bod- 
ies of the two galaxies coalesce, but are still discrete entities, which 
also provides us with the relaxed requirements in terms of spatial 
resolution (see Fig.[T]l. In addition, the relatively low surface bright- 
ness extended gas produced during such events (see, e.g., the tidal 
tail in Fig. [T) will also allow us to explore the ability of MOAO 
in recovering such faint discrete features. Trying to detect these 
features again emphasizes our choice of simulations and is really 
leading toward constraints that optimize for surface brightness de- 
tection and away from resolving structures within disks and in the 
chaotic central regions of advanced mergers. In the next three sec- 
tions, we describe the different steps of this study: (a) Simulation 
of data-cubes (Sect. 3); (b) Simulations of MOAO PSFs (Sect. 4); 
(c) Kinematic analysis of the data-cubes (Sect. 5). 



3 AN END-TO-END SIMULATION PIPELINE FOR 3D 
SPECTROSCOPY 

3.1 General Description 

The main steps of the process can be summarized as follows. First, 
a data-cube with with the spatial resolution of the telescope diffrac- 
tion limit (i.e., ~ A/2D, where D is the telescope diameter) is gen- 
erated. For each pixel of this high resolution data-cube, a spectrum 
is constructed from observations of local galaxies or results from 
numerical simulations. In the second step the spatial resolution of 
the data-cube is reduced by convolving each spectral and spatial 
pixel of the high resolution data-cube by a PSF. This PSF used in 
this convolution is representative of the optical path through the at- 
mosphere and the telescope up to the output of the (optional) AO 
system. In the third step, the spatial sampling of the data-cube is 
reduced to match that of the IFU of the simulated instrument. Fi- 
nally, realistic sky as well as photon and detector noise are added. 
Note that this software can also be used to simulate images of dis- 
tant galaxies, since an image can be viewed as a particular case of a 
data-cube, with a single broad-band spectral channel. We now de- 
scribe each step in detail; all input parameters are summarized in 
Tabled 

3.2 High resolution data-cube 

At very high spatial resolution, emission lines with kinematics 
driven only b y gravitational motions is well described by a sim- 
ple Gaussian teeauvais & Bothuij[l999l) . Under this assumption, 
for each pixel of the data-cube, only four parameters are required 
to fully define a spectral line. The first three parameters are the 
position in wavelength, the width, and the area (or, equivalently 
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Figure 1. Morphology and kinematics of the two objects used in tlie simulations. Up: UGC5253 (rotating disk). Down: Merging between two Sbc galaxies. 
From left to right: continuum map, monochromatic map, velocity field, and velocity dispersion map. The original spatial scale is ~85 pc/pix for UGC5253, 
and ~200 pc/pix for the merger simulation. These templates are rescaled in terms of size and flux to simulate realistic distant galaxies. 



the height) of the emission line. The current version of the soft- 
ware only models emission lines and does not take into account 
the detailed shape of the continuum in galaxies - it is simply mod- 
eled as a constant in f,\. So only one parameter is required to set 
the level of this pseudo-continuum around the emission line. We 
assumed a perfect Atmospheric Dispersion Corrector (ADC) and 
did not take into accou nt atmospheric refraction effects (see, e.g., 
lOoncharov et al.ir2007l) . Each spectrum is generated in the obser- 
vational frame, at a given spectral sampling of A„i,J2R, where R 
is the spectroscopic power of resolution of the instrument, and 
'lobs ~ (1 -I- z)Aem, -^em being the rest-frame wavelength of the emis- 
sion line, and z the redshift of the simulated object. During this pro- 
cess, the rest-frame line width is multiplied by (1-l-z), as one needs 
to take into account the widening of emission lines with increasing 
redshift. 

All four parameters can be extracted from observations of lo- 
cal galaxies. In the following, we use Fabry-Perot (FP) observa- 
tions of the Ha emission dist ributions of nearby galaxies obtained 
as part of the GHASP survey dAmram et al.l2003) . From these data, 
we can extract four parameters, wavelength, width, area, and con- 
tinuum level to construct the velocity field, the velocity dispersion 
map, the flux map of the Ha emitting gas, and the continuum map 
of the galaxy. The software first re-scales all these maps at a given 
angular size (in arcsec) provided by the user. These maps are then 
interpolated at a spatial sampling of ~ A/2D. This spatial sampling 
is motivated by the fact that the AO PSFs used have been simulated 
at this sampling (see §4). The software then re-scales the overall 
amplitude of the continuum map at a given integrated number of 
photons using an integrated magnitude hiab directly provided by 
the user. This magnitude is converted to the number of photons 
per spectral pixel depending on the telescope diameter D, the in- 
tegration time ti„,g, and the global transmission /,„„., of the system 
(atmosphere excluded). The overall amplitude of the Ha map is 
also rescaled with a given integrated number of photons, derived 
from the integrated continuum value and a rest-frame equivalent 
width EWq provided by the user, the latter being re-scaled in the 
observed-frame by multiplying by (1+z). 

All parameters can also be extract ed from outputs of hydro- 
dynamical simulations jCox et aL I l2004h . In this case, the last two 
parameters (area and continuum level) can respectively be extracted 
from total gas and stellar surface density maps which are by- 



products of the numerical simulation, also rescaled in terms of size 
and flux. 

3.3 Modeling the IFU and the detector 

Each monochromatic slice of the high resolution data-cube is con- 
volved by a PSF with matching spatial sampling. This PSF must be 
representative of all elements along the optical path, from the at- 
mosphere to the output of a n (optional) Adaptive Optic system. Be- 
cause the isoplanetic patch ( lFriedll98ll : the median value at Paranal 
is ~2.4 arcsec atA~ 0.5 //m, which leads, e.g., to ~10 arcsec a.lA~ 
1.6 fim) is larger than the individual FoV of the IFU (typically a 
few arcsec, depending on the size of objects at a given redshift), the 
same PSF can be used to convolve the data-cube regardless of posi- 
tion within the IFU. We also neglected the variation of the PSF with 
the wavelength (with a FWHM varying as in a Kolmogo rov 
model of the atmospheric turbulence, see, e.g.. lRoddieilll98ll) , as 
we are only interested in the narrow spectral range around a single 
emission line. 

The next step is to reduce the spatial sampling of the data- 
cube. This is done by re-binning each monochromatic channel of 
the data-cube at the pixel size A,„i of the simulated IFU. A wave- 
length dependent atmospheric absorption curve taken from ESO 
ParanaQ is then applied to each spectrum of the data-cube. Sky 
continuum, detector dark level and bias are then added to the spec- 
tra. We used a sky spectrum model (including zodiacal emission, 
thermal emission from the atmosphere, and an average amount of 
moonlight) from Mauna Ke^ which has the advantage over other 
available sky spectra to be very well sampled with 0.2 A/pixel. Pho- 
ton and detector noise are then added to each individual exposure. 
The detector noise is due, in the NIR, to the dark current (dark), and 
readout noise (ron). The simulation pipeline generates ndit data- 
cubes with individual exposure time of dit, which are combined by 
estimating the median of each pixel to simulate several individual 
realistic exposures. Since we have only included random noise, it 
is similar to having dithered all of individual exposures and com- 
bining them after aligning them spatially and spectrally. The spec- 
troscopic SNR is then derived as follows: 

' www.eso.org/observing/etc 

2 www.gemini.edu/sciops/ObsProcess.obsConstraints/ocSkyBackground.html 
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0{i^,jy,k,i) = 



^0{ix, jy, k^) + S (( J, jy, + ron- + dark 



where 0(i_^,jy,k,i) and S(i_^,jy,A) are respectively the object and 
sky flux per dit (after accounting for atmospheric transmission) in 
the spatial position (i.v, /,) of the data-cube (in pixels), and at the 
spectral position k,^ along the wavelength axis (in pixels). In the fol- 
lowing, the "maximal SNR in the emission line in the pixel (i^, jy)" 
refers to MAX^t^ [SNR{ix, jy, k^)], and the "total SNR" refers to the 
flux weighted average of this quantity over the galaxy (i.e., the spa- 
tially integrated SNR in the spectral pixel where the emission line 
peaks, and not the spatially integrated SNR within the whole emis- 
sion line width, as sometimes used in other studies). 



4 GENERATING MOAO PSFS 

The coupling between the MOAO system and the 3D spectrograph 
is captured through the MOAO system PSF. Therefore, it is a cru- 
cial element that needs to be carefully simulated, and cannot be ap- 
proximated by, e.g., a simple Gaussian. We describe here how these 
PSFs were modelled, and briefly describe how the MOAO correc- 
tion influences the PSF shape at the spatial scale of the instrument 
IFU. 



4.1 MOAO PSFs for EAGLE 

A preliminary study of the PSF shape generated by a hypothetical 
MOAO system can be found in Neichel et al. (2006). In the fol- 
lowing, we used eight PSFs generated in a similar manner as in 
the Neichel et al. study. Briefly, three off axis guide stars, located 
at the edges of an equilateral triangle, are used to perform a to- 
mographic measurement of the turbulent atmospheric volume. We 
have used a turbulent profile typical of that on C erro Paranal, whic h 
can be modeled using three equivalent layers jFusco et alj|l999h . 
The seeing and the outer scale of the turbulence were respectively 
set to 0.95 arcsec (at 0.5 fim, the ~0.75 percentile on Paranal) and to 
22m (the median value on Paranal). Guide stars are assumed to be 
bright (V~13), i.e., WFS noise is neglected. The optimal correction 
is deduced from the characteristics of the turbulence volume and 
applied assuming a single Deformable Mirror (DM) per direction 
of interest, here taken as the center of the guide star constellation. 
To explore a wide range of correction, we consider different GS- 
constellation sizes, as well as a DM with an inter-actuator size of 
either 1 m or 0.75 m in the pupil plane (for a 42 meter telescope), 
which corresponds to 42x42 or 56x56 actuato rs (see Tab. [3). Of 
note, the SPHERE eXtreme A O (XAO) project jSeuzit et ahHoO^ 
will use a 41x41 actuators DM IFusco et alj2006h . PSFs are shown 
in Figure |2] 

Usually, AO PSFs are characterized using the Strehl Ratio 
[SR], which provides a useful way for comparing different PSFs 
relatively to the diffraction limited case. However, MOAO does 
not provide diffraction limited PSFs but performs only partial cor- 
rections in order to increase the spec troscopic coupling with the 
3D spectrograph jAssemat et alj|2007h . which cannot be easily de- 
scribed by using the SR. Instead, the spectroscopic coupling can be 
directly quantified by the fraction of light under the PSF that en- 
ters a spatial element of resolution of the IFU, i.e., the Ensquared 
Energy entering an aperture equal to twice the IFU pixel size. Such 
a parameter is an extrapolation of the classical "Encircled Energy" 
used in optics to characterize optical system quality, and is directly 



Pitch 


FoY WFS 


EE 


Strehl Ratio (%) 


1.00 


4.00 


22 


0.1 


1.00 


3.00 


24 


0.2 


1.00 


2.00 


26 


0.6 


1.00 


1.00 


31 


3.9 


1.00 


0.50 


34 


6.7 


1.00 


0.25 


40 


11.5 


1.00 


0.00 


43 


10.5 


0.75 


0.00 


47 


11.7 



Table 3. AO system parameters used for the simulation of the EAGLE 
PSFs. The EE is given in 0.15 arcsec, and the pitch corresponds to the 
inter-actuator distance on the telescope pupil, in meters. The second col- 
umn gives the diameter of the constellation of three stars (at equal distance 
of the central galaxy) used to sense the wavefront, in arcmin. Strehl Ratios 
are given in the last column. 



related to the achieved SNR. It is therefore a more natural choice 
to parametrize the MOAO performance than the SR. In the pres- 
ence of speckle noise on the PSF as it is the case in the end-to-end 
MOAO simulations used here, the SR measurement can be sev- 
erally affected by this noise, which produces artificially low SR 
values and contributes to disconnect the SR measurement from the 
spectroscopic coupling: at a given EE, a PSF affected by speckle 
noise will have a lower SR than a PSF without speckle noise, all 
else being equal. We list both SR and EE for the simulated PSFs in 
Tab. [3] 

To explore how the AO correction impacts the PSF shape, we 
first d erived the PSF FWHM using Sextractor (Berlin & Arnout^ 
1 19961) . We compared these to their Ensquared Energy [EE] in a 
aperture equal to twice the pixel size (i.e., 0.15 arcsec) in Fig. [3] 
The EE is measured as the fraction of light under the PSF that en- 
ters this aperture. Fig.|3]reveals two regimes, depending on the EE. 
In the first regime, the EE and the FWHM are roughly linearly 
correlated, so that the EE can alternatively be used as an indirect 
measure of the FWHM. This regime corresponds to the formation 
of a central diffraction limited core in the PSF. When the adaptive 
system does not correct for the overall atmospheric turbulence, it 
is indeed well known that the resulting PSF can be approximately 
described by a double core-halo structure, with a coherent central 
diffraction-limited peak surrounded by seeing-limited wings. This 
regime corresponds to an improvement of the spatial resolution, as 
long as the PSF FWHM remains larger than twice the IFU pixel 
scale. When the PSF FWHM becomes smaller than twice the IFU 
pixel size, the spatial resolution is now set up by the IFU pixel scale. 
In this case, the AO corrections does not dramatically influence the 
spatial resolution any longer, but can still provide better EE and 
SR by further narrowing the core below the IFU pixel scale (see 
Tab. O, which translates into an improvement of the SNR. When 
EE> 27-30%, the diffraction limited core is now well formed, which 
explains why the measurement of the FWHM of the PSF saturates, 
since the FWHM is then only sensitive to the PSF central core as 
it does not "see" the residual halo around it. In this new regime, 
the SR keeps rising with the EE because the AO system is still able 
to bring energy located in the PSF halo into the central diffraction 
limited peak. In most situations, the MOAO system will provide 
EE larger than this limit, which justifies the choice of characteriz- 
ing the EE in an aperture equal to twice the pixel size. 

In the intermediate regime where the PSF FWHM is smaller 
than twice the IFU pixel scale, but where the diffraction core is not 
completely formed (e.g., EE ranging between 23 and 35% in Fig. 
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PaTtiiTictcr 


Description (Unit) 


EAGLE/E-ELT 


FALCONA'LT 


Ml 


Telescope Primaiy min'or size (m) 


42 


8.2 


M2 


Telescope Secondary mirror size (m) 





1.116 


R 


Spectral resolution power 


5000 


5000 




IFU pixel size (arcsec) 


0.075" 


0.125" 


hraiism 


Telescope + Instrument transmission factor 


0.2 


0.2 


dark 


Dark level (e/sec/pix) 


0.01 


0.01 


ron 


Read-out-noise (e/pix) 


2.3 


2.3 


dit 


detector integration time (s) 


3600 


3600 


ndit 


Number of dit per simulation 


200-24-8 


200-24-8 


z 


Redshift of the object 


4.0 


1.6 


S 


Object angular diameter (arcsec) 


0.8 


2.0 


mAB 


Integrated magnitude of the object continuum (AB mag) 


24.5 


22.5 




Emission line wavelength (A) 


3727 


6563 


EWo 


Equivalent width at rest of the emission line (A) 


30 


50 



Table 2. First and second columns: Input parameters of the simulation pipeline provided by the user Third column: parameters used for the simulations of 
z=4 galaxies with EAGLE on the E-ELT. Fourth column: parameters used for the simulations of z= 1 .6 galaxies with FALCON on the VLT. 



the PSF halo can be an important limitation in the observations 
of high redshift galaxies, because it can result in a mixing of the 
spectra coming from adjacent spatial elements of resolution. In the 
following, we use the total energy entering the first ring of 0.15 
arcsec apertures around the one centered on the PSF to measure 
this effect. This provides a useful proxy to quantify the amount of 
light coming from adjacent spatial element of resolution, and pol- 
luting the element of resolution centered on the PSF. We called this 
effect "PSF cross-talk" (PCT), by analogy with the instrumental 
crosstalk, which measures the amount of scattering that spreads in- 
formation beyond the PSF due to instrumental effects. Figure[3]also 
compares the PCT with the EE for the set of eight simulated PSFs, 
and shows a linear correlation between these two parameters: while 
the AO system provides a better correction, the energy located in 
the wings of the PSF is b rought into the central diffraction-limited 
core (see lAssemat et al1 [2007). increasing the EE and decreasing 
the PCT. PCT gives as useful measurement of the "PSF contrast": 
having a PCT smaller than 50% means that the light entering a 
spatial element of resolution is not dominated by polluting light 
coming from adjacent spatial elements of resolution. In the EA- 
GLE case and at a pixel scale of 75 mas, this occurs at EE~35% 
(see Fig.[3]l- 



4.2 MOAO for FALCON 

An extensi ve study of the FAL CON MOAO system has been car- 
ried out by lAssemat et alj ( 12007) . As part of this study, a simula- 
tion pipeline of FALCON PSF was developed. We used this tool to 
generate a set of eight PSFs (see Fig.|4l(. All PSFs were simulated 
using typical atmospheric conditions at ESO/Paranal. Briefly, a me- 
dian seeing of 0.81 was used, with an outer scale of the turbulence 
of 24m, and an atmosphere model composed of three layers located 
at [0,l,10]km with respective weights of [20, 60, 20]%. We set the 
parameters of the MOAO system in order to explore a wide range 
of EE in 0.25 arcsec (twice the pixel size), ranging from 19% to 
46% (s ee Tab.[4ll. Technica l details about the concept can be found 
in, e.g., |Puech et ai]( l2005h . 

For this set of PSFs, we compared their FWHM and PCT to 
their EE in 0.25 arcsec in Figure[5] which reveals a similar behavior 
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Figure 2. Transverse cut of the simulated EAGLE MOAO PSFs. The EE 
in 0.15 arcsec increases from bottom to top: 22%, 24%, 26%, 31%, 34%, 
40%, 43%, and 47%. The vertical dashed-lines represent the spatial element 
of resolution used for EAGLE simulations, i.e., 0.15 arcsec. 



compared to Fig. [3] Note that when EE>30%, the FWHM becomes 
lower than 0.25 arcsec (i.e., twice the pixel size). This means that 
when EE>30%, the spatial resolution is oversampled at twice the 
IFU pixel size. This is why, in the following, all EE are computed 
in a 0.25 arcsec size aperture. Finally, PCT becomes smaller than 
50% for EE~35%, which is the limit for a good "PSF contrast" in 
the FALCON case, at a pixel scale of 0.25 arcsec. 
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Figure 3. Top: FWHM vs. EE (in 0.15 arcsec) for the eight simulated EA- 
GLE PSFs. The dotted Hne represents the diffraction limit on a 42m ELT 
in H band, i.e., Adiff(H) = 0.008 arcsec. The black line is an exponential 
fit between the FWHM and the EE, and the dotted-dashed line is a simple 
linear fit to the three first points. The dashed Unes shows when the I^WHM 
becomes lower than 0.15 arcsec, i.e., twice the pixel size (under-sampling 
of the PSF). Bottom: PCT vs. EE for the eight sunulated EAGLE PSFs. The 
black line is a simple linear fit between PCT et EE. 



Table 4. AO system parameters used for the simulation of the FALCON 
PSFs. The EE is given in 0.25 arcsec, and the value quantifying the correc- 
tion corresponds to the radial number n of Zemike modes corrected. The 
third column gives the diameter of the constellation of three stars (at equal 
distance of the central galaxy) used to sense the wavefront, in arcsec. Fol- 
lowing Puech et al. (2005), n=7 corresponds roughly to a ~1.4m pitch, and 
n=9 to a ~l.lm pitch. The first PSF is the seeing limited case, with no 
correction. Strehl Ratio are given in the last column. 



Correction FoV^fs EE Strehl Ratio (%) 






N.A. 


19 


1.03 


7 


300 


26 


1.56 


7 


180 


27 


1.51 


7 


60 


31 


2.26 


7 


30 


34 


4.04 


7 


10 


35 


4.55 


7 


60 


39 


20.82 


9 


60 


46 


30.03 



Figure 4. Transverse cut of the FALCON MOAO PSFs simulated. The EE 
in 0.25 arcsec increases from bottom to top: 19%, 26%, 27%, 31%), 34%, 
34%, 35%, 39%, and 46%. The vertical dashed-Hnes represent the spatial 
element of resolution used for FALCON simulations (0.25 arcsec aperture) 
within which EE are derived. 
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5 SIMULATING AND MEASURING 2D KINEMATICS 
5.1 Kinematic measurements 

The simulated data-cubes are analyzed using a data analysis 
pipeline similar to those generally used to analyze data of high red- 
shift galaxies. During this process, each spatial pixel of the simu- 
lated data-cube is fitted with a Gaussian in wavelength, whose po- 
sition and width correspond to the velocity and velocity dispersion 
of the gas in this spatial pixel. Because the accuracy of these mea- 
surements is driven by the flux within the emission line relative to 
the noise in the continuum, one has to define a kinematic signal-to- 
noise ratio SNRkin specific for this measurement. We chose to rely 



Figure 5. Top: FWHM vs. EE for the eight simulated FALCON PSFs. The 
left most point corresponds to the seeing limited PSF. The horizontal line 
represents the diffraction limit on the VLT in H band, i.e., Xdiff(H) = 0.04 
arcsec. The dashed lines show when the spatial resolution, as measured by 
the IJWETM, is lower than 0.25 arcsec, i.e., twice the pixel size. The black 
line is a linear fit between IJWHM and EE (the last point, where the I^WHM 
saturates, has been discarded from the fit). Bottom: PCT vs. EE for the eight 
simulated FALCON PSFs. The left most point corresponds to the seeing 
limited PSF. The black fine is a linear fit between the PCT and the EE. 
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on the definition of lFlores et alj l l2006h . who defined S NRi^j,, as the 
total flux in the emission line divided by the noise on the continuum 
Ccmiunuum and V^,„.i . the number of pixels within the emission line. 
We emphasize that this SNRkm is diff'erent from the classical SNR, 
as defined in Sect. 3.3: the latter quantifies the signal to noise ratio 
relative to the sky and detector noise (detectability of the emission 
line), while the former characterizes the capability of measuring 
the first and second order mo ments of th e emission line, relative to 
the noise in the continuum ( Sarzi et al.ll2006) . Although they are 
different in nature, we find that both SNRs roughly correlates pro- 
vided that SNR > 3 (corresponding the SNRm,, ~3-4). Because a 
lot of data-cubes are simulated and need to be analyzed, one has to 
rely on a fully automatic analysis pipeline. Tests show that pixels 
with a "kinematic" SNRki,, < 3 have very large errors mostly due to 
misidentification of emission lines with noise peaks. Hence, only 
pixels having a 5 NRi,i„ of at least three have been kept for further 
analysis. 

5.2 Kinematic diagnostic diagrams 

To help distinguish between the two kinematic templates (i.e., the 
major merger and the rotating disk), we fitted a rotating disk model 
to each velocity field. This model assumes a constant dynamical 
center, PA, inclination, and Kns o\er the galaxy, and a simple arc- 
tan shape for the rotation curve. During the fit, the pixels of the 
velocity field are weighted by the corresponding error in bins of 
SNR, which were determined using Monte-Carlo simulations. Fi- 
nally, the residuals between the velocity field and the fitted model 
are derived, together with the reduced map. 

We then adopted two criteria for judging whether or not the fit 
to a particular simulated galaxy is appropriate. Because we are in- 
terested only in large-scale motions (i.e., circular rotation vs. non- 
structured motions), we need to adopt a criterion not affected by 
smaller-scale variations due to non-circular perturbations like, e.g., 
spiral arms or bars. For characterizing the 2D-fit as a whole, we 
therefore rely on the median ^,^,(VF) over the galaxy. When smaller 
than unity, it means than at least 50% of the pixels on the velocity 
field are reliably fitted by a rotating disk model, and that galaxy 
motions are actually dominated by circular rotation. In addition, 
we choose to compare the distance between the fitted dynamical 
center and the velocity dispersion barycenter. We used such a crite- 
rion because in the case of data which are heavily smoothed by the 
seeing, a process that leads to more complex dynamics like merg- 
ing would only rarely have its dynamical center of the entire system 
coincident with the maximum in the distribution of velocity disper- 
sions. For this, we derived the associated uncertainty, using classi- 
cal propagation of errors due to the uncertainties on the dynamical 
center fit and on the sigma barycenter. We then derived a param- 
eter AcTpfnj., defined as the distance between the fitted dynamical 
center and the velocity dispersion barycenter, normalized by the 
corresponding uncertainty. Having Ao"p,„^-j < 1 means that the fit- 
ted dynamical center is found close enough to the sigma barycenter 
location relative to the associated uncertainty, to be considered as 
compatible with the dynamics of a rotating disk. 

In summary, the kinematics of a galaxy can be considered as 
compatible with that of a rotating disk, provided that it falls on a 
region of a "diagnostic diagram" (xl,(y^),^crpeai!s) where the me- 
dian normalized residual is smaller than one, and provided that 
peaks ^ 1- We need to emphasize that this classification (as any) 
is not without limitations. As an example, let us consider a merg- 
ing pair of galaxies. If both progenitors are equally distant from 
the center of mass of the system, the velocity dispersion map will 



show two peaks equally distant from this center, and the resulting 
barycenter of the sigma map may then be found very close to the 
dynamical center of orbital motions between the two progenitors if 
the separation of the two is small compared to the intrinsic resolu- 
tion of the data and the exact distribution of the emission line gas. 
Of course, that is a limitation of most techniques as it can be diffi- 
cult to identify close merging pairs at high redshift. It is important 
to be aware that there are some configurations that are difficult to 
distinguish between. Be that as it may, we use this approach as a 
first attempt in quantifying the rotation in very distant galaxies. 



6 SIMULATIONS 
6.1 Assumptions 

It was decided not to simulate observations in the K-band as the 
detectability of lines in the thermal infrared is highly dependent 
on the achieved emissivity of the instrument. To generalize these 
first simulations and to remove additional free parameters, such as 
the telescope and instrumental emissivities and temperatures, we 
consider the H band only, assuming no thermal emission from the 
instrument and the telescope in the H-band. Our tests show that 
this assumption leads to overestimate the SNR in the H-band by 
no more than ~10%. This makes these first simulations less depen- 
dent of the telescope design (e.g., number of mirrors), environmen- 
tal conditions (site selection), and instrument characteristics (e.g., 
number of warm mirrors). Moreover, working in the H-band will 
help disentangle spatial effects from flux limitation arising from 
the high thermal background in the K-band. 

For this study, we assumed a spectral power of resolution 
i?=5000 as a compromise between our desire to minimize the im- 
pact of the OH sky lines and not wanting to over-resolve the line 
by a large factor. Moreover, appropriate targets in the NIR are usu- 
ally selected as they have emission lines that fall in regions free 
of strong OH lines. We estimated that about one third of the H- 
band is free from strong sky background variations (i.e., larger 
that 10% of the sky continuum) in continuous windows of at least 
200 km/ill Therefore, for simplicity, all simulations have been per- 
formed without OH sky lines, accounting for the sky continuum 
only. Usually, R=3000 is considered as a strict minimal spectral 
resolution to work between OH sky lines. In comparison, R=5000 
will allow a better sky subtraction and limit kinematic measurement 
(i.e., velocity and velocity dispersion) uncertainties. In addition, for 
simplicity, we assume that the [Oil] emission line is a single line, 
instead of a doublet. This does not influence any results presented 
here, as these are scaled with the total flux, i.e. the flux inside both 
lines of the doublet (assuming the low density limit for the ratio of 
the two lines). 

The instrumental parameters of the simulation were set to typ- 
ical values used in the IR (see Tab|2j, relying on a cooled Rock- 
well HAWAII-2RG IR array working at ~80K, as described by 
iFinger et al. !(|2003). An array of this type is already used in SIN- 
FONI, and arrays with similar characteristics will be implemented 
in next VUT generation IR instruments such as HAWK-I, KMOS, 



^ Tills obviously depends on the spectral resolution and sampling of the 
sky background. These figures were determined using the Mauna Kea sky 
model Introduced In Sect. 3.3, which has a spectral resolution of 0.4 A and 
Is sampled at the Nyquist rate. This gives a spectral power of resolution 
R~22500 at ^=9000A. 
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or X-SHOOTER. Table[2]summarizes the main parameters used in 
the simulations. 

Previous studies of FALCON have considered a pixel scale of 
0.125 arcsec tPuech & Sa yede 2004; Assemat et al. 2007). We will 
thus adopt this scale in the following. For EAGLE, we chose to 
firstly explore a pixel scale of 75 mas. This choice is suggested by 
recent 2D-kinematical studies of distant galaxies with GIRAFFE, 
which show that a relatively coarse pixel scale can already provide 
us with enough information to properly recover large scale mo- 
tions in distant galaxies, and then d iscriminate between simple ro- 
tating disks and major mergers (see lFlores et alJ[2006l : |Puech et al.l 
l2006ah . 



6.2 Simulation of z~4 galaxies with EAGLE 

We conducted simulations at z=4, where the [Oil] emission line 
falls in the H band. Both objects have been rescaled in half-light 
radius an d flux to account for e volu tion and redshif t (see bottom 
of Tabllll. iFerguson etal] ( l2004h and iBouwens et al] | |2004|) found 
that the typical r;,a// of z~4 galaxies is ~0.2 arcsec in the UV. 
Hence, we adopted an optical diameter of 0.8 arcsec for z~4 galax- 
ies ( Rppt ~ Irkaif). The continuum magnitude was chosen consid- 
ering [^shMaeT^ 12006), who studied z~4 Lyman Break Galax- 
ies luminosity functions in the UV. They found a typical UV mag- 
nitude of 24.5. We neglected any color corrections between UV 
and blue optical ([OII]/i3727A): as spiral Spectral Energy Distri- 
butions [SEDs] rise from UV to optical, this choice can be consid- 
ered as a minimal value for a typical pseudo-continuum for z~4 
galaxies. We assumed an integrated rest-frame equivalent-width 
EWo([OII] )=30A, extrapolating the median values found in z<l 
galaxies bv lHammer et al.l ( ll997b . This leads to an integrated emis- 
sion line flux ~7.5xI0""*erg/s/cm^ (observed-frame). 



6.3 Simulations of z~1.6 galaxies witli FALCON 

We choose as a rest-frame Equivalent Width EWo for the Hg emis- 
sion li ne, the minimal value observed in the sample of lErb et"an 
1 I2OO6I) at z~2. In their sample, the continuum at 1500A ranges 
from m^B = 22 to 26. Assuming a U -H color of 3.5, typical of Sa- 
Sb galaxies (Mann ucci et al. II2OO1I), and assuming that this value 
can be extrapolated at z~1.6, this leads us to a minimal magni- 
tude in the H band of ~22.5. As a cross-check, w e looked at the 
sample of BM/BX galaxies of lReddvetZI ( I2OO6I) . We found 22 
sources having a redshift between 1.4 and 1.8, for which Reddy 
et al. give I and K bands AB magnitudes. Assuming that the H- 
band magnitude can roughly be approximated by an average be- 
tween J and K bands magnitudes, we found in these 22 sources, a 
median (mean) H-band magnitude of 22.55 (22.56). We can thus 
assume with some confidence that mAB(H)=22.5 is quite repre- 
sentative of luminous z~1.6 galaxies. This leads to an integrated 
emission line flux ~4.9x 10~'^er g/s/cm^ (observed-frame). Finally, 
following [Ferguson et a i] l l2004h . we assumed a typical half-light 
radius r;,„/y of 0.5 arcsec i n UV (~ 1 500 A) for z~1.6 galaxies. A 
similar value was found bv lDahlen e t al. ( 2007) at 2800A, with 0.4 
arcsec for z= 1.75 galaxies withM, ^IO'^Mq. We took 2 arcsec as a 
linear diameter for z~1.6 galaxies, neglecting the k-morphological 
correction between UV and rest-frame optical (Ha). The parame- 
ters adopted for the simulations are listed in Tab.|2] 



7 RESULTS 

7.1 Tlie Design Requirements for EAGLE on tlie E-ELT 

As explained above, the effects of the PSF on the data-cube are 
twofold. On one hand, the FWHM (and the PCT) influences the 
spatial resolution of the observations. On the other hand, the EE 
strongly influences the number of photons reaching the detector. 
Because the EE roughly correlates with the FWHM (and the PCT), 
both effects are related. An easy way to distinguish them and to 
focus on spatial features is to simulate data-cubes with very large 
SNRs. We therefore first ran a set of simulations with t„„g ~ 200h: 
Figure [7] shows the velocity fields, velocity dispersion maps, and 
emission line flux map for the z~4 simulated rotating disk. The ef- 
fect of the PSF can clearly be seen on each of these three maps 
(from left to right). More flux (SNR) if brought toward the galaxy 
center, the "spider" shape of the isovelocities are more robustly re- 
covered, and the peak at the dynamical center of the dispersion map 
is not so heavily smoothed as the correction increases. These effects 
are even clearer in the simulation of the major merger: the two 
peaks in the emission flux maps (see also the velocity dispersion 
maps), which correspond to the two progenitors, are clearly better 
distinguished as the correction improved, as are the non-circular 
motions in the velocity field. The two progenitors can be distin- 
guished provided that EE > 26%. It is worth to note that this EE 
corresponds to a FWHM smaller than 2 pixels (see Fig. O, i.e. to 
a regime where the central diffraction peak is already formed (see 
Sect. 4.1). As the two progenitors are separated by only 2-3 pixels, 
a poor correction mixes the two corresponding peaks in the emis- 
sion line flux maps and/or the velocity dispersion maps. Obviously, 
detecting smaller spatial-scale features would require a finer pixel 
scale. However, at a given pixel size, these simulations suggest that, 
from a purely spatial point of view (i.e. at infinite SNR), an optimal 
way to couple AO with 3D spectroscopy is to tune the MOAO sys- 
tem in order to get a PSF whose FWHM is not larger than twice the 
pixel size. Improving the correction further will not provide a better 
spatial resolution, but at that point a finer pixel scale may help (see 
Sect. 4). However, even with a relatively coarse sampling, a higher 
EE, will in turn lead to a higher SNR in a given integration time: 
with realistic integration times, the optimal choice of AO correc- 
tion must take into account achieving sufficient SNR to detect low 
surface brightness regions. 

To investigate the impact of the SNR on the AO requirements, 
we ran simulations with shorter total integration times: we adopted 
perhaps more realistic values of 24 and 8 hours. With an integration 
time of 24hr, the total SNR over the entire galaxy ranges between 
~11 and ~22, depending on the EE (see Fig. llOt. This SNR goes 
down between ~7 and ~ 12 with 8hr of exposure time, scaling as the 
square root of the integration time, as expected. For a given surface 
brightness distribution and integration time, the SNR scales linearly 
with EE, as expected in a background dominated regime. Adopting 
a total SNR of 5 as a lower limit to recover reliable and meaningful 
information over a galaxy, one can derive that between 1.5 and 4 
hours of integration time is a strict minimum to recover z=4 galaxy 
kinematics, depending on the EE provided by the MOAO system 
and the pixel scale. 

In Fig. [8] we plot radial-average profiles of the maximal SNR 
reached in the emission line for the rotating disk simulations, both 
with 24 and 8hr of exposure time, and with 22 and 47% of the EE 
(the smallest and largest values used in the simulations). This figure 
shows the effect of increasing the EE for a given integration time: 
it mainly increases the SNR in the central part of the galaxy (which 
has some structure fine enough to benefit from good AO correc- 
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Figure 6. Kinematics of a major merger as derived from simulations of EA- 
GLE observations at z~4, witli 75mas/pixel, t„„g=24h and EE=34%. The 
velocity fields (first column), the velocity dispersion maps (second column), 
and the emission line flux maps (last column, units are in median counts per 
dit) are showfn. The object size is 0.8 arcsec in diameter, which represents 
-0.12 kpc atz=4. 

tion). On the other hand, the outer and lower surface brightness 
regions are much more difficult to recover: with 8hr of exposure 
time, emission lines are only detected over only a limited region of 
the galaxy; at least 24hr are required to obtain full coverage of the 
galaxy or merger, i.e., a SNR of at least 3 in the outer parts. We find 
that at least 24hr and an EE of at least 34% are needed to recover 
the merger's low surface brightness tail (see Fig. [Sjl. We note that 
this EE corresponds to the EE required to have a PCX smaller than 
50%, i.e., to have a good "PSF contrast", which limits the spectro- 
scopic coupling between adjacent spectra (see Fig.O. 

Even when the full kinematics of the galaxy or merger is not 
recovered, the total SNR is always larger that 5, down to 8hr of 
exposure time. We see that this partial information is sufficient to 
recover the dynamical nature of both the major merger and the ro- 
tating disk, for a large range of EE (see Fig.|9}. Although based on a 
limited number of cases, this suggests that it is not mandatory to re- 
cover the full 2D Idnematics to infer the dynamical nature of distant 
galaxies. The dynamical state of a galaxy is mostly reflected in its 
larg e scale motions, ra ther than in its smalf scafe perturbations (see 
also lFlores et al ]|200 6). This is why it appears sufficient to recover 
only part of the 2D information, providing that this information is 
reliable enough (i.e., recovered with enough SNR). We will discuss 
this point further in Sect. 8.1. 

7.2 The Design Requirements for FALCON on the VLT 

We proceed as described in §3.6, and ran a set of simulations with 
a total integration time t,„,j=200h (see Fig. lilt . To distinguish a 
merger from a simple rotating disk in both the cr and velocity maps, 
we find that a minimum ~30% of EE (in 0.25 arcsec) is required. 
Again, this EE corresponds to a FWHM smaller than 2 pixels (see 
Fig.H}. 

With 24hr of exposure time, the total SNR ranges between ~6 
and ~11, and between ~4 and ~8 with 8hr, depending on the EE 
(see Fig.ll4t. Reaching a total SNR of 5 thus requires at least 8hr 
with a minimal EE of 30%, or more than 3hr, assuming an EE of 
46%. 

We find that recovering the full kinematic information requires 
at least 24hr of exposure time and an EE of 35% (see Fig. 112b . 
Again, this EE corresponds to a PCT smaller than 50 % (see Fig.[ 5ll, 
i.e., to a good "PSF contrast". In a previous study, lAssemat et al.l 
l l2007l) found that 35% of EE was required for detecting the Hor 
emission line at a 3(T level. However, one can expect that for mea- 
suring higher moments of the emission line, i.e., the first and second 
moments, a larger SNR will be required. If the velocity measure- 
ment does not seem to be affected much by the SNR, the veloc- 
ity dispersion clearly suffers from a too low SNR, since the two 
progenitors can hardly be directly identified in Fig [12] relying on 
velocity dispersion maps only. 
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Radius [pix] 

Figure 8. Radius-average profile of the maximal SNR in the [Oil] emission 
line of the rotating disk, obtained after 24hr (full lines) and 8hr (dotted 
lines) of exposure time. In each case, the upper curve is obtained for the 
highest EE (47%), and the bottom one for the smallest EE (22%). Error- 
bars account for the 1-sigma dispersion in SNR and positions within each 
ring of the profile. Only pixels having SNRkin >3 have been considered to 
build the profile. 
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Figure 10. Flux- weighted SNR in the center of the emission line, integrated 
over the galaxy for EAGLE simulations with 75mas/pix and 24 or 8 hr of 
exposure times (respectively: full lines and dotted lines). The black curves 
correspond to the rotating disk, while the red ones correspond to the major 
merger. 
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Figure 7. Kinematics of a rotating disk (tliree top lines) and a major merger (tliree lowest lines) as derived from simulations of EAGLE observations at z~4, 
with 75 mas pixel"' and t,„,j=200h. The velocity fields (VF), the velocity dispersion maps (o"), and the emission line flux maps ([Oil], units are in median 
counts per dit) are shown. Each column corresponds to a given EE in 0.15 arcsec, as indicated at the top. 






Figure 9. From left to right: Diagnostic diagram for identifying rotating disks in EAGLE simulations of z~4 galaxies with 75mas/pixel and t„„g ~200h, 24h, 
and 8h. Each point represents a different AO correction, from the worst one (EE=22%, see the biggest point) to the best one (EE=47%, see the bigger square). 
The black track corresponds to the rotating disk, and the red one to the major merger. 



With 8hr of integration time, it is possible to recover the true 
underlying nature of dynamics of high redshift galaxies, using the 
diagnostic diagram (see Fig. |13t . However, a minimal EE of 26% 
and 35% are respectively necessary for integration times of 24 and 
8h. However, the diagnostic diagram generally fails when the total 



SNR is smaller than ~5 (see Fig.ll4t. confirming once again that 5 
appears to be the lowest acceptable SNR. 
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Figure 11. Kinematics of a rotating disk (tliree top lines) and a major merge (three lowest lines) as derived from simulations of FALCON observations at 
z~1.6, with t,„,g=200h. For both objects, the velocity fields (VF), velocity dispersion maps (cr), and emission line maps (Ha, units are in median counts per 
dit) are shown. Each column coiTesponds to a given EE in 0.25 arcsec, as indicated at the top. 



8 DISCUSSION & SUMMARY 

8.1 Spatial "scale-coupling": Influence of the pixel size for 
EAGLE 

We also ran simulations of the EAGLE instrument concept for the 
E-ELT with a smaller pixel scale, 50 mas (see Fig. |16t . FWHMs 
and PCTs are shown in Fig. [15] The two progenitors of the ma- 
jor merger can be clearly separated in both the cr and the [Oil] 
maps (see Fig.|16t with EE >15% (in 0.1 arcsec), which again cor- 
responds to the EE where the FWHM becomes smaller than two 
pixels. 

As expected, the achieved total SNR at a given integration 
time and EE, is lower with a 50 mas pixel scale than with one of 



75 mas pixel"' (see Fig.|18t. Reaching a minimal total SNR of 5 
requires at least 8hr of exposure time with an EE of 21%, or at least 
~3hr with a larger EE. 



With 50 mas pixel"', we find that the low surface brightness 
tidal tail of the merger can be detected in t;„,j,=24h with at least 
32% of EE. In contrast to the previous cases, this EE does not cor- 
respond to the one giving 50% of PCX: this condition now gives 
a much smaller EE (smaller than 10, see Fig. 1151), which does not 
provide enough SNR to make this distinction (Fig.llSt. With the 75 
mas pixel" ' scale, detecting the tidal tail required a total SNR of 
~15: this is roughly the total SNR that is now provided by 32% of 
EE with the 50 mas pixel"' scale. This is obviously related to subdi- 
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Figure 13. Diagnostic diagram for identifying rotating disks in FALCON simulations of z~1.6 galaxies with 125mas/pixel and t,„,g ~200h, 24h, and 8h (from 
left to right). Each point represents a different AO correction, from the worst one (EE=19%, see the biggest point) to the best one (EE=46%, see the bigger 
square). The black track corresponds to the rotating disk, and the red one to the major merger. 
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Figure 16. Kinematics of a rotating disk (thi'ee top lines) and a major merger (thi'ee lowest lines) as derived from simulations of EAGLE observations at z~4, 
with 50mas/pixel and t,„,g=200h. The velocity fields (VF), the velocity dispersion maps (cr), and the emission line flux maps ([Oil], units are in median counts 
per dit) are shown. Each column corresponds to a given EE in 0. 1 arcsec, as indicated at the top. 



viding the total flux of the tail into smaller and smaller sampling of 
a structure with a relatively smooth surface brightness distribution. 

Using the diagnostic diagram, it is possible to disentangle the 
major merger from the rotating disk with a pixel scale of 50 mas 



pixel ' (see Fig. I17t . Again, the diagnostic diagram tends to fail 
when the total SNR is smaller than 5 (see Fig.|18ll. 

The comparison of the 50 and 75 mas pixel"' scales reveals 
that the former does not bring significant more information com- 
pared to the latter, relatively to the goal of distinguishing between 
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Figure 17. From left to right: Diagnostic diagram for identifying rotating disks in EAGLE simulations of z~4 galaxies with 50 mas pixel"' and t„„j ~200h, 
24h, and 8h. Each point represents a different AO correction, from the worst one (EE=12%, see the biggest point) to the best one (EE=37%, see the bigger 
square). The black track corresponds to the rotating disk, and the red one to the major merger. 
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Figure 12. Kinematics of a rotating disk and major merger as derived from 
simulations of FALCON observations at z~1.6 (0.125 arcsec/pixel), with 
t„„j=24h and EE=35%. For both objects, the velocity fields (first column), 
velocity dispersion maps (second column), and emission line maps (third 
column, units are in median counts per dit) are shown. The objects size is 2 
arcsec in diameter, which represents ~0.24 kpc at z=1.6. 



the major merger and the rotating disk. This is related to the fact 
that, as mentioned before, the dynamical state of a galaxy is re- 
flected in its large scale motions, rather than in its small scale per- 
turbations. Hence, having a very fine pixel scale does not provide 
any advantage in this respect, and might limit the SNR achievable 
for a given integration time. Obviously, such a conclusion depends 
on the spatial scale of the kinematic feature to be recovered: if 
one wants to recover smaller-scale kinematic details like, e.g., the 
precise rotation curve, or to detect clumps in distant disks, a finer 
spatial scale will be then preferable. This optimal "scale-coupling" 
between the IFU pixel scale and the spatial scale of the physical 
feature that one wants to recover using this IFU can be quantified 
by the ratio between the size of this feature (here, the galaxy di- 
ameter, as one wants to recover the large-scale rotation) and the 
size of the IFU resolution element. At z~4, this "scale-coupling" 
is found to be 5.3 using the 75 mas pixer' scale, and 8, using the 
50 mas pixer' scale. In other words, our simulations suggest that a 
"scale-coupling" of 5.3 is already large enough to recover large- 
scale motions in z~4 galaxies. This is in line with 3D observa- 
tions of z~0.6 galaxies with FLAMES/GIRAFFE, which provide s 
a "scale-coupling" of about 3 dFlores et alj200dlYang et al.l2008h : 
this is the minimum value necessary to ensure that each side of the 
galaxy is at least spatially sampled by the IFU at the Nyquist rate. 
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Figure 14. Total flux-weighted SNR over the galaxy for FALCON simu- 
lations with 24 and 8 hr of exposure times (respectively: full lines, dashed 
lines, and dotted lines). The black curves correspond to the rotating disk, 
while the red ones correspond to the major merger. 



Using this minimum ratio, one finds that a pixel scale of at least 
~130 mas is mandatory for studying large-scale motions in z~4 
galaxies if they are simply scaled versions of nearby galaxies. On 
the other hand, the smallest pixel scale that can be used for a given 
scientific goal is limited by the achievable SNR at this scale, which 
is in turn constrained by the "PSF contrast" as discussed above in 
Sect. 7. 



8.2 Influence of Telescope Diameter 

As discussed previously, in most cases MOAO will provide 3D 
spectroscopy of distant galaxies with a spatial resolution limited 
by the IFU pixel scale, because MOAO essentially provides partial 
AO corrections. Everything else being equal, using a larger aper- 
ture telescope will then only influence the integration time needed 
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Figure 15. Top: FWHM vs. EE for the eight simulated EAGLE PSFs. The 
dotted line represents the diffraction limit on a 42m ELT in the H-band, 
i.e., Ajiff(H) = 0.008 arcsec. The dashed hnes shows when the FWHM 
becomes lower than 0.1 arcsec, i.e., twice the pixel size (under-sampling). 
Bottom: PCT vs. EE for the eight simulated EAGLE PSFs. The black line 
is a linear fit between FWHM et EE. 
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Figure 18. Total flux-weighted SNR over the galaxy for EAGLE simula- 
tions with SOmas.pix"' and 24 and 8 hi' of exposure times (respectively: 
full lines, dashed lines, and dotted lines). The black curves con'espond to 
the rotating disk, while the red ones correspond to the major merger 



to reach a given SNR, thanks to the larger collecting power of the 
primary mirror. 

We show in Sect. 4, that an MOAO-fed IFS on the VLT could 
provide us with improved 2D-kinematics of z<2 galaxies. In Fig. 
II9I we show what the maps of the same galaxy would look like 
if observed with EAGLE on the E-ELT, using 50 mas pixel"' and 
t;„,g=2h. Observations, with 24% of EE in 0. 1 arcsec, provide a ra- 
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Figure 19. Kinematics of a rotating disk (first line) and a major merger (sec- 
ond line) as derived from EAGLE simulations at z~1.6, with t,„,g=2h and 
EE=24% in 0.1 arcsec. For both objects, the velocity fields (first column), 
velocity dispersion maps (second column) and emission line maps (third 
column) are shown. Object size is 2 arcsec in diameter, which represents 
~0.24 arcsec at z=1.6. 



dial SNR profile comparable to the one obtained with FALCON 
after 24hr of exposure time and EE=39% (see Fig. I20t . The E- 
ELT, equipped with an multi-object IFS, even with a low EE, would 
represent a huge improvement over 8-meter class telescopes, even 
equipped with future-generation MOAO-fed instruments. In return, 
the VLT, equipped with an MOAO-fed spectrograph would pro- 
vide us with an important advancement in our knowledge of dis- 
tant galaxies, increasing the number of spatial element of resolu- 
tion from typically 2-3 to 13-20. This will be possible at the VLT 
with the Laser Guide Star Facility and SINFONI (pixel scale of 50 
mas). This is already availabl e with OSIRIS at Keck (with 20 to 10 
mas pix AO-corrected; see iLaw et alj|2007l : IWright et al.ll2007h . 
Other VLT second-generation NIR spectrographs, such as KMOS, 
will on the other hand allow us to gather further information on 
kinematics of 24 distant objects simultaneously, which is a great 
advantage in obtaining the statistical properties. However, it only 
has seeing limited performance over a 0.2 arcsec pixel"' scale. High 
spatial resolution integral field spectroscopy at the VLT will remain 
limited to single object observations. The number density of ob- 
jects at z ~ 1.4-2.5 is expecte d to be ~ 9 per square arcmin, down to 
Iab -25 ( ISteideletal.1l2004t) . KMOS will undoubtedly identify in- 
teresting objects amount the ensemble of za:2 galaxies with which 
to observe in more detail with AO-fed integral field spectrographs 
on 8-m class telescopes. However, an MOAO-fed spectrograph on 
currently available or a future extremely large telescope would op- 
timally take advantage of these high target densities. 



8.3 Summary: High Level Design Requirements for MOAO 
systems 

We first presented simulations of z~4 galaxies with very large 
SNRs, to focus on spatial smoothing effects due to partial correc- 
tions of the atmospheric turbulence by the MOAO system. In all 
cases, we find that the EE required to clearly identify the two com- 
ponents of the major merger considered in our simulations roughly 
corresponds the EE that provides us with a FWHM smaller than 2 
pixels. For z=4 galaxies and using a 75mas pixel scale, this corre- 
sponds to an EE of ~25% in 150 mas, while this condition translates 
into 15% in 100 mas with a 50mas pixeL' scale. At z=1.6, an EE of 
at least ~30% within 0.25 arcsec is required, using a 0.125 arcsec 
pixel"' scale. Considering the spatial resolution only, (i.e., with an 
infinite SNR), these specifications are indications of the minimum 
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Figure 20. Radius-average profile of the maximal SNR in the emission line, 
obtained after 24hr of exposure time and 39% of EE in 0.25 arcsec with 
FALCON (full lines), and after 8hr of exposure time and 24% of EE in 0. 1 
arcsec with EAGLE (dotted lines). Error-bars account for bootstrap uncer- 
tainties on the medians in SNR and positions within each ring of the profile. 



we will realistically require to recover spatially resolved kinematic 
information on distant galaxies. 

Considering the impact of signal to noise, if one wants to re- 
cover the full 2D kinematics of z=4 galaxies, at least 24 hr of in- 
tegration time, and 34% of EE (in 0.15 arcsec, using the 75 mas 
pixer' scale) are needed (or 32% in 0.1 arcsec using the 50 mas 
pixel"' scale). For z=1.6 FALCON simulations, at least 24 hr and 
35% of EE in 0.25 arcsec are needed for recovering the full 2D 
kinematics in our simulations. When the SNR provided by the 
MOAO system is large enough, these EE roughly corresponds to 
a PCX smaller than 50%. This provides a good "PSF constrast" at 
a given spatial scale, which avoids significant amount of polluting 
light from adjacent spectra. 

At lower SNR, the required EE naturally depends on the sur- 
face brightness of the faintest regions to be detected, as well as on 
the adopted pixel scale. The choice of the latter strongly depends 
on the scientific goal and the size of the features to be recovered by 
3D spectroscopy: a good "scale-coupling" between them must pro- 
vide a least a few spatial element of resolutions per important scale. 
Because the dynamical state of galaxies is mostly reflected in their 
large-scale motions, this means that only a few spatial element of 
resolution are needed per galaxy diameter: we performed additional 
simulations (not shown here) using another rotating disk with a less 
steep rotation curve (UGC6778, see lGarrido et al.l2002l) , and found 
no quantitative difference in our results using this other rotating 
disk. Although based on a limited number of simulations, we find 
that indeed a spatial resolution element of 150 mas (pixel size of 75 
mas) seems to be a good choice for recovering large scale motions, 
providing relatively large SNR with a "scale-coupling" ~5.3 ele- 
ment of resolution per galaxy diameter (0.8 arcsec). At the 50mas 
pixer' scale, a 42 m ELT will typically need 8 hr of exposure time 
to reach a SNR of at least 5 over a galaxy. Such a SNR allows us 
to use the diagnostic diagram as a useful tool to recover the dy- 
namical nature of rotating disks and major mergers. Because of the 



much smaller collecting area of the VLT, this exercise is limited 
to 8 hr with at least 35% of EE in 0.25 arcsec, for z=1.6 galaxies. 
In this case, the "coupling factor" is 4. These results are consistent 
with previous studies at lower redshift (z~0.6), where the dynam- 
ical nature of galaxies have been recovered using GIRAFFE, with 
only 3 or so resolution elements (at 0.52 arcsec pixel"') across each 
galaxy. 



9 CONCLUSION 

We have developed software capable of end-to-end simulations of 
integral field spectrograph fed by an MOAO system. We have used 
this software to investigate a limited number of simulations in H- 
band, in order to give insights into how such MOAO-fed systems 
should perform on the VLT and the E-ELT. By requiring that the in- 
strument is able to distinguish between a rotating disk and a major 
merger, we were able to constrain the required Ensquared Energy 
to be able to recover their large-scale motions. Separating the pro- 
genitors of the major merger requires only modest EE (15-26% at 
z=4 in a sample of 100-150mas; 30% at z=1.6 in 0.25 arcsec), but 
long exposure times (~24hr). Higher EEs are generally needed to 
recover the full 2D kinematics in our simulations, typically 35%. 
Provided that the total SNR is large enough over the galaxy (typ- 
ically ~ 5), it is possible to use more sophisticated methods, such 
as the diagnostic diagram introduced here, to distinguish between 
both types of objects, even at relatively low integration time for 
such distant objects (8 hrs). To generalize these results beyond what 
we described here, we have related these specifications to the con- 
cept of "PSF contrast", which is the amount of polluting light from 
adjacent spectra and which drives the minimal EE at a given spa- 
tial scale. The choice of the latter is driven mostly by the "scale- 
coupling" which is the relative size of the IFU pixel and the size of 
the features that the observe wishes to recovered in the data cube. 
Distinguishing between a major merger and a disk, for example, 
can be largely done by investigating the large scale motions within 
a system, and a relatively coarse spatial resolution appears then to 
be sufficient. Given this situation, a pixel scale of 50-75 mas seems 
to be a relatively good choice, since it provides sufficient SNR in 
less time than a finer sampling and also has the advantage of relax- 
ing the MOAO system requirements. More stringent requirements 
can be set by attempting to resolve and investigate structure in high 
redshift galaxies and that will be simulated in subsequent papers. 
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